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ABSTRACT 

We have computed the gravitational wave signal from supernova core collapse using the presently 
most realistic input physics available. We start from state-of-the-art progenitor models of rotating 
and non-rotating massive stars, and simulate the dynamics of their core collapse by integrating the 
equations of axisymmetric hydrodynamics together with the Boltzmann equation for the neutrino 
transport including an elaborate description of neutrino interactions, and a realistic equation of state. 
Using the Einstein quadrupole formula we compute the quadrupole wave amplitudes, the Fourier wave 
spectra, the amount of energy radiated in form of gravitational waves, and the signal-to-noise ratios 
for the LIGO I and the tuned Advanced LIGO ( "LIGO II" ) interferometers resulting both from non- 
radial mass motion and anisotropic neutrino emission. The simulations demonstrate that the dominant 
contribution to the gravitational-wave signal is produced by neutrino-driven convection behind the 
supernova shock. For stellar cores rotating at the extreme of current stellar evolution predictions, the 
core-bounce signal is detectable (S/N > 7) with LIGO II for a supernova up to a distance of ~ 5 kpc, 
whereas the signal from post-shock convection is observable (S/N > 7) with LIGO II up to a distance 
of ~ 100 kpc, and with LIGO I to a distance of ~ 5 kpc. If the core is non-rotating its gravitational 
wave emission can be measured with LIGO II up to a distance of ~ 15 kpc (S/N > 8), while the 
signal from the Ledoux convection in the deleptonizing, nascent neutron star can be detected up to a 
distance of ~ 10 kpc (S/N > 8). Both kinds of signals are generically produced by convection in any 
core collapse supernova. 

Subject headings: supernovae: general — hydrodynamics — relativity — gravitational waves 



1. INTRODUCTION 

For more than two decades astrophysicists struggle to 
compute the gravitational wave signal produced by core 
collapse supernovae (for a review, see eg., Miiller 1997). 
If the core collapse and/or the resulting supernova 
explosion involves time-dependent asphericities such 
that the third time-derivative of the quadrupole moment 
of the mass-energy distribution is nonzero, part of the 
gravitational binding energy liberated in the event 
will be emitted in the form of gravitational waves. 
Such non-sphericities can be caused on large scales 
by the effects of rotation JMiilled" Il982t iFinn fc Evans! 
19901 lM5nchmever ef.aZl 119911 lYamada fc Satol I199E " 



Zwerger fe Mulled 119971 IRampp. Miiller fc Ruffert 



1998; IDimmehneier. Font fc Mull 

" Frver. Holz fc Hughes! I2002t IFrver et 

Imamura. Pickett fc Durisenl l2003t lOtt et.al\ | 200£ 



and low-mode convection ijHerantl 1995; 
iScheck et. al 2003), on small scal es by flow fluctuations 
due to high-mo de convection ijMuller fc Jankal Il997t 
iFryer e£.al||20p 3|). and by anisotropic neutrino emission 
(lEpsteinlll97ft iBurrows fc Haves! Il99"fl iMiiller fc Jankal 
1997). 

Theoretical predictions of the gravitational wave 
signal from core collapse supernovae have been, and still 
are, hampered by the complex nature of the supernova 
explosion physics which is not yet fully understood 
(jBuras et.a l 2003) . Reliable simulations require realistic 



(rotating) pre-collapse stellar models, the incorporation 
of a realistic equation of state (EoS), a detailed modeling 
of weak interaction processes, Boltzmann neutrino trans- 
port, multi-dimensional hydrodynamics, and relativistic 
gravity. However, all past studies aimed at studying the 
gravitational wave signature of core collapse supernovae 
have considered greatly simplified parameterized models 
involving one, several or all of the following approxi- 
mations: A polytropic equation of state, a simplified 
description of weak interactions and neutrino transport 
or none at all, parameteriz ed pre-collaps e stellar models, 
and Newtonian gravity llMulled Il982t iF inn fc Evan; 
19901 iMonchmever etoflj 19911 lYamada fc Satol 1199 
" Zwerger fc Mulled 119971 IRampp. Miiller fc Rufferti 
19981 iDimmelmeier. Font fc Mulled 120021 

Frver. Holz fc Hughes! 12002} IFrver et.al\ 



Kota ke. Yamada fc 
2003|) . In addition. 



Imamura. Pickett fc Durisenl |2003; 

Satol 120031 IQtt et.al\\2(M IShibatal 

most of these simulations focussed 
on the effects of rapid rotation, and on the bounce signal 
covering only the evolution up to a few ten milliseconds 
after core bounce. 

From observations as well as theoretical modeling it 
is now commonly accepted that core collapse super- 
novae do generically involve asphericities besides those 



expected in the case of a progenitor with a sufncientl 
rapidly rotating core (for a 

(1998)). These asphericities are impor tant or may 



rapidl y rotating core (for a review see, eg., [Miiller 
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| 1994f iBurrows. Haves fc Frvxelll 119951: Uanka fc MuTIer! 
1996;Hn^^SH2003). Recent state-of-the-art axisym- 
metric (2D) si mulations using a Boltzmann solver for 
the v transport l)Rampp fc J anka 2002; Bu ras et.al 2003) 
and simplified (as far as neutrino transport and neutrino- 
matte r interactions are con cerned) three-dimensional 
(3D) l)Frver fc Warrenl l2002j) simulations both confirm 
that convective overturn indeed occurs in the v heat- 
ing region and is helpful for shock revival, thus mak- 
ing explosions p ossible even when sph erically symmet- 
ric models fail ijJanka fc Mtillerl 11996ft . However, de- 
spite of strong convective action and a corresponding en- 
hancement of the efficiency of neutrino energy transfer 
to the post-shock matter, even the up to now most real- 
istic simulations of both non-rotating and r otating pro- 
genit or models do not produce explosions ijBuras et.al\ 
2003). Multi-dimensional simulations further demon- 
strate that vigorous aspherical motion is unavoidable in 
the slowly deleptonizing and cooling proto-neutron star 
llKeil. Janka fc Mullerlll99fifl Both the hot bubble con- 
vection ijMuller fc Jankall997[) and the convection in the 
proto-neutron star (see below) as well as the resulting 
anisotropic neutrino emission ijBurrows fc Haves! 11 996: 
iMuller fc JankalH997j) give rise to an interestingly large 
gravitational wave signal. 

Concerning consistent pre-collapse stellar models 
from evolutionary calculations of rotating stars, ma- 
jor progress has recently b een a chieved by the work of 
iHeger. Langer fc Wooslevl i|2000l) . who have performed 
the up to now most realistic evolutionary calculations 
of rotating massive stars to the onset of core collapse 
including the effects of rotation, mixing, transport of an- 
gular momentum, most recently, also of magnetic torques 
l)Heger eta l 2003). Including the latter the central an- 
gular velocities become of the order of O.lrad/s only. 
Despite of these fairly slow rotation rates the cores are 
still too fast to lead to typically expected neutron star 
natal spins, and angul ar momentum loss is considered 
as a serious problem l|Wooslev fc Hegerl EOQSh The low 
iron core rotation rates, however, suggest that in spite 
of angular momentum conservation during core collapse, 
triaxial rotation-triggered instabilities need not to be ex- 
pected in the collapse simulations. For the time being, 
we therefore constrain ourselves to axisymmetric models. 

In the following we present the gravitational wave sig- 
natures of the currently most elaborate and detailed core 
collapse supernova simulations. One of the considered 
models starts from an iron core with a rotation rate at 
the onset of collapse that is in the ballpark of predic- 
tions from the latest generation of stellar evolution mod- 
els. The simulation also contains a state-of-the art de- 
scription of neutrino-matter interactions and Boltzmann 
neutrino transport, a realistic equation of state, and 2D 
axisymmetric hydrodynamics taking into account gen- 
eral relativistic effects. We also present the gravitational 
wave signature of a state-of-the-art non-rotating super- 
nova model. This model, which is based on another, 
less massive progenitor, serves the purpose to demon- 
strate that also from non-rotating stars interesting grav- 
itational wave signals can be expected, which are strong 
enough to be detectable for supernovae that occur even 
at large distances in our Galaxy. Finally, we present 
the gravitational wave signal expected from the Ledoux 
convection in a proto-neutron star at times after the ex- 



plosion has taken off (~ Is after core bounce). 

The paper is organized as follows. In the next section 
we discuss the physics employed in our models, and the 
numerical techniques used to perform the simulations. 
In Section 3 we present the method used to extract the 
gravitational wave signature, and in Section 4 we describe 
the results of our investigation. Finally, in Section 5 we 
discuss the implications of our findings, and give some 
conclusions. 

2. INPUT PHYSICS AND NUMERICAL TECHNIQUES. 

For integrating the equations of hydrodynamics we em- 
ploy the Newtonian finite-volum e code PROMETHEUS 
i|Frvxell. Miiller fc Arnettl 11989ft . This second-order, 
time-explicit Godunov code is a direct Eulerian implc- 
mentation of the Piecewise Parabolic Method (PPM) 
l)Colella fc Woodward! Il984). and is based on an exact 
Riemann solver. 

We use spherical polar coordinates (r,9,ip), and as- 
sume axial symmetry. However, motions in the in- 
direction are allowed, and the azimuthal (rotational) ve- 
locity may be nonzero. In the rotating model the rota- 
tion axis coincides with the symmetry axis at 9 = 0. In 
this model we have also imposed equatorial symmetry. 
The 2D computational grid consists of 400 logarithmi- 
cally spaced zones in radial direction and, if equatorial 
symmetry is assumed, of 64 equidistant angular zones 
covering an angular range < 9 < tt /2. If equatorial 
symmetry is not assumed, the angular grid consists of 
128 equidistant zones in the angular range < 9 < n. 

The calculat ions are performed with th e equation of 
state (EoS) of lLattimer fc Swestvl (|1991f) using an m- 
compressibility of bulk nuclear matter of 180 MeV. This 
EoS is based on a compressible liquid-drop model includ- 
ing nuclei, nucleons, e~, e + , and photons. 

Neutrino transport is computed wit h the Boltzmann 
solver sche me describe d in d etail in iRampp fc Janka! 
( 2002) and iBuras et.al\ l|2003fl for v and v of all three 
flavors. For each angular bin of the numerical grid the 
monochromatic moment equations for the radial trans- 
port of v number, energy, and momentum are solved. 
This set of equations is closed by a variable Eddington 
factor that is calculated from the solution of the Boltz- 
mann equation on an angularly averaged stellar back- 
ground. Beyond this "ray-by-ray" approach we also take 
into account the coupling of neighboring ra ys by lateral 
advec tion terms and v pressure gradients ijBuras et.al\ 
2003). For the transport an energy grid of 17 geomet- 
rically spaced bins is used with centers from 2 MeV to 
333 MeV. 

The simulations were done with a Poisson solver that 
calculates the gravitational potential from the 2D mass 
distribution. General relativistic effects are treated ap- 
proximately by modifying the spherical part of the grav- 
itational potential with correction terms due to pres- 
sure and energy of the stellar medium and neutrinos, 
which are deduced from a comparison of the Newto- 
nian and r elativistic equations of motion in spherical 
symm etry l)Keil. Janka fc Mullerlll996T: IRampp fc Janka! 
2002). The v transport contains gravitational redshift 
and time dilation, but ignores the distinction between 
coordinate radius and proper radius. Comparison with 
fully relativistic, one-dimensional simulations showed 
that these approximations work well at least when the 
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deviat ions of the metric coeffic ients from unity are mod- 
erate l|Liebendorfer etaZll2003l). 

We considered two (solar metallicity) progenitor mod- 
els with mai n sequence masses of 11.2 M Pl a nd 15 Mq, 
respectively l)Wooslev. Heger fc Weaverl l2b02"l. The for- 
mer progenitor was adopted as given by the stellar evo- 
lution calculation (model sllnrl80), while angular mo- 
mentum was added to the 15 Mq progenitor (model sl5r; 
see also Buras et.al2003). The angular frequency of 
model sl5r was assumed to be constant in the Fe and 
Si core and to decrease like r -3 / 2 outside of 1750 km 
(1.43 Mq). This rotation profile is in agre ement with 
predi ctions from stellar evolution models ijHeger et.al\ 
2003). We have adopted a value of 0.5 sec -1 for the cen- 
tral angular frequency, which corresponds to an initial 
value of (3 = 0.05%, where (3 is the ratio of the rota- 
tional energy and the (absolute) gravitational binding 
energy. This is a relatively small rotation rate compared 
to those typically assumed in parameter studies based o n 
rotating polytropes (see, eg., iZwerger fc Mulled ITool : 
iDimmelmeier Font fc Mullerl (j2002j) : lOtt et.aiTfiOO^ ). 
However, it is of the order of the largest rotation rates 
predicted by state-of-the-art evolutionary calculations of 
massive stars, and even too fast with respect t o obse rved 
rotation rates of young pulsars t'H eger etai\ l)2003l) . see 
also Sect. EJ- Hence, as far as the state of rotation is 
concerned, our model is meant to represent not an ex- 
treme, rare case, but is in accord with the conditions in 
stellar iron cores as expected from current models for the 
evolution of rotating massive stars. Due to the smallncss 
of the assumed rotation rate we regarded it unnecessary 
to relax the initial model into rotational equilibrium (the 
ratio between the centrifugal force and the gravitational 
force is below 1% in the progenitor model). 

We have also analyzed a proto-neutron star model 
(pns!80), which is the 180° analogue of the model dis- 
cussed in lKeil. .Tanka fc Miillerl 1|199(tD . which covered an 
angular region of only 90° size. Model pnsl80 provides 
complementary information about the gravitational wave 
signal expected from core collapse supernovae, namely 
the wave emission from the long-lasting (> 1 sec) non- 
radial motion driven by the deleptonization of the form- 
ing neutron star. The ~ 1.1 M© convective core of 
the proto-neutron star of model pns!80 was formed in 
the c ore collapse of a 15 M Q star (iKeil. Janka fe M iillcr 
11996ft . Its evolution was followed by iKeill l)1997[) for 
a period of more than 1.2 sec with a variant of the 
PROMETHEUS hydrodynamics code (see above) assum- 
ing axial symmetry and including a flux-limited (equi- 
librium) neutrino diffusion scheme, which was applied to 
each angular bin separately ("1.5D"). The computational 
grid consisted of 100 radial zones extending out to an ini- 
tial radius of 58.8 km and 120 equidistant angular zones 
covering an angular range < 9 < n, ie., no equatorial 
symmetry was assumed. During the simulated evolution 
the grid was contracted with the neutron star shrinking 
to a final radius of about 20 km l|Keil. Janka fc M iillcr 

mm. 

3. CALCULATION OF THE GRAVITATIONAL WAVE 
SIGNATURE 

The (quadrupole) gravitational wave amplitudes and 
energy spectra resulting from a nisotropic mass motion 
are computed as described in Mi iller fe Jankal l|1997fl 



(eqs. [10]-[12]) using the Einstein quadrupole formula 
in the numerically convenient form ulation derived by 
iBlanchet. Damour fe Schaferl (^990), and by standard 
FFT techniques. The wave amplitudes due to the 
anisotropic neutr ino emission are o btaine d with the for- 
malism given in iMiiller fc Jankal 119971) fen s. [28]- [31]) 
which is based on the work of lEpsteinl 1)1978). 

In order to find out whether the gravitational wave sig- 
nals of our model stars (assumed to be located at some 
given distance) are detectable by terrestrial gravitational 
wave experiments we have also computed the correspond- 
ing (optimal with respect to the relative orientation of 
source and detector) signal-to-noise ratios S/N obtained 
from matched filtering, which are determined by the re- 
lation ijFlanagan fe HughesllT998|) 

(W = 4f&, (!) 

where h{y) is the Fourier transform of the gravitational 
wave amplitude h{t) and Sh(v) [Hz -1 ] is the spectral den- 
sity of the strain noise in the detector. Eq. JQ) can also 
be written as 



/>oo 


' \h{y)\ v 1 / 2 ' 


/o 





dlogv (2) 



where /i rms (i/) [Hz~ 1/< 2 1 = \/ Spjv) is the spectral strain 
sensitivity (see, eg., iFritschefl ((2002)), which gives the 
rms noise level of the antenna at the frequency v. 
The second expression for the signal-to-noise ratio S/N 
(Eq.EJ is useful when discussing the spectra of the grav- 
itational wave signals of our models (see Section 4.2). 

In particular, we have determined the signal-to- 
noise ratios of our models both for the LIGOI and 
the tuned Advanced LIGO ("LI GOH") interfe romet- 
ric gravitational wave detectors ijFritschell 12002]) . For 
LIGOI we used a fit to the strain sensitivity curve 
LIGOI SRD Goal available from the LIGO web site 
www.ligo.caltech.edu/docs/G/G030014-00/. In case of 
the tuned Advanced LIGO interferometer we have com- 
puted signal-to-noise ratios for the three different strain 
sensitivity curves shown in Fig.^ 

Concerning the sensitivity curves shown in this plot 
and some of our other figures, we point out that both 
the curves for LIGO II as well as the curve for LIGO I 
are for a single interferometer. In the case of LIGOI, 
there are two nominally identical 4 km instruments, plus 
a 2 km instrument at Hanford. This increases the net 
sensitivity by a factor \J~2 in the ideal case of statistical 
independence and gaussian noise. For LIGO II, there are 
3 interferometers (two at Hanford, one at Livingston). 
They could all be tuned the same way for a directed 
search, or tuned differently. If all are tuned the same way, 
the net sensitivity is a factor better at all frequencies 
than that of a single interferometer. The signal-to-noise 
ratios given in this publication are for the complete de- 
tector consisting of two interferometers (LIGO I) or three 
like-tuned interferometers (LIGO II), respectively. 

Assuming axial symmetry and an observer located 
at an angle with respect to the symmetry axis of 
the source, the dimensionless gravitational wave ampli- 
tude h(t) is related to the quadrupole wave amplitude 
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Fig. 1. — Strain sensitivity curves of the tuned Advanced LIGO 
interferometer ("LIGO II") used in our analysis. The solid curve 
has an increased sensitivity at ~ 1 kHz and a rather wide trough 
of near-maximal sensitivity, while the dotted curve has a peak sen- 
sitivity of 1.2 10 — 24 at about 700 Hz, an d ali gns with the peak in 
energy in some of our models (see Section l4.2l . The third (dashed) 
sensitivity curve is tuned for sources which emit more at low fre- 
quencies (the so-called NS-NS-tuned instrument). 



(measured in units of cm), the lowest-order non- 
vanishing term of a multipolc expansion of the radia- 
tion field into pure-spin tensor harmonics (see Eq. [9] of 
Miiller 1997), according to 



1 



I A E2 
R 



(3) 



where R is the distance to the source. In the following 
we always assume sin 2 6 = 1. 

4. RESULTS 
4.1. Wave amplitudes 

The quadrupole wave amplitudes of models sl5r, 
sllnrl80 and pnsl80 are shown in Figs.|2| 01 and ^ re- 
spectively. Besides the waveform resulting from the con- 
vective flow, each figure shows the signal due to the 
anisotropic neutrino emission. Inserts further illustrate 
details of the wave signal for chosen periods of time, and 
demonstrate that the wave signal consists of many quasi- 
periodic variations on time scales of a few milliseconds 
which remain unresolved when the whole time evolution 
of the models is shown. 

The quadrupole wave amplitude Af 2 (Eq.EJ) of model 
sl5r at bounce is of the order of 20 cm showing (see 
insert of Fig.l2"|l the typical prominent spike with sub- 
sequent ring-down of a type I gravitational wave signal 
of a rotating core, which bounces due to the stiffening 
of the equation of state around nuclear matter density 
ijMonchmever eta71ll991l) . The post-bounce wave am- 
plitude stays below ~ 10 cm until about 150 msec when 
convective mass motions driven by neutrino heating in 
combination with the increasing rotation rate of the post- 
shock region amplify the wave signal. The growing rate 
of rotation of the post-shock region is caused by the ac- 
cretion of matter with higher angular momentum (as the 
initial rotation profile is constant to the edge of the Fe 
core) through the shock, and the specific angular mo- 
mentum increasing even beyond the edge of the iron 
core. The post-shock region starts to exhibit large-scale 
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Fig. 2. — Gravitational wave quadrupole amplitude Ajn 2 vs. time 
(post bounce) due to convective mass flow and anisotropic neutrino 
emission (thin line) for the rotating delayed explosion model sl5r of 
IBuras et.al 12003). The insert shows an enlargement of the signal 
around the time of bounce. 
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Fig. 3. — Same a s Figure l2l but for the delayed explosion model 
sllnrlSO of Buras et.al (2003). The insert shows an enlargement 
of the signal between 150 msec and 175 msec after core bounce. 



non-radial pulsations (Fig.[5j) which produce a long-time 
(~ 10 — 20 msec) variability on top of the rapidly varying 
(^rnsec) wave signal of up to 130 cm amplitude. 

The non-rotating, and hence initially spherically sym- 
metric model sllnrl80 shows, as expected, no signal at 
bounce (Fig. 01 • Until 50 msec after bounce the wave am- 
plitude of this model remains tiny (< 1cm). Only when 
the neutrino-driven convection sets in at around 40 msec 
to 50 msec, the wave amplitude grows to the size of a 
few centimeters. Eventually, as the convective activity 
becomes stronger, the maximum (absolute) amplitudes 
rise to values of several 10 centimeters. 

Both models, the rotating (sl5r) as well as the non- 
rotating one (sllnrl80), do not explode within the dis- 
played time intervals. After an explosion has started we 
expect a slow decay of the convective activity behind the 
expanding supernova shock, and hence a declining gravi- 
tational wave amplitude from this region over a period of 
a few 100 msec. The gravitational wave signal afterwards 
will be dominated by the emission from convection inside 
the cooling nascent neutron star. 

The proto-neutron star model exhibits a gravitational 
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Fig. 4. — Same as Fig,l2l but for the proto-neutron star model 
pnslSO of IKeill 119971) . Note that no bounce occurs in this model, 
which follows the neutrino cooling phase of the nascent neutron 
star. Time is measured from the start of the simulation at about 
20 msec after stellar core collapse. 




Fig. 5. — Four snapshots of the entr opy distributio n of the ro- 
tating delayed explosion model sl5r of Buras et.al (2003). The 
side length of the plots is 600 km, and the numbers in the top left 
corners give post-bounce times. Bright regions are high entropy 
bubbles, the shock is visible as the sharp, deformed discontinu- 
ity, and the flattened proto-neutron star is the low-entropy oblate 
ellipse at the center. 



wave signal consisting of rapid (~ msec) quasi-periodic 
temporal variations of roughly constant amplitude (~ 
cm) modulated on time scales of ~ 100 ms (Fig.QJ). While 
the former variations reflect mass motions, and non- 
radial proto-neutron star oscillations occurring on the 
hydrodynamic (sound crossing) time scale, the long-time 
signal modulations are caused by the varying strength 
of the activity due to Ledoux convection in the proto- 
neutron star (FigJSJ. 

The wave amplitude associated with anisotropic neu- 
trino emission shows much less time structure and slower 
temporal variation than that produced by non-radial 
mass motion. It is characterized by an overall rise of 
the wave amplitude in case of the rotating model sl5r, a 
quasi-periodic oscillation with a time scale of ~ 50 msec 




5 10 15 20 
R [km] 

Fig. 6. — Relative variation of the proton-to-baryon ratio (elec- 
tron number fraction) due t o convectiv e mass flow in the proto- 
neutron star model pnslSO of Kcil ( 1997) 0.7sec after core bounce. 
Proton-rich matter rises, while neutron-rich matter sinks inward. 



in case of the non-rotating model sllnrl80, and an even 
slower temporal variation (0.2 — 0.4 sec) in case of the 
proto-neutron star model pnsl80. 

In case of the rotating model (sl5r) the temporal vari- 
ation of the quadrupole wave amplitude of the neutri- 
nos reflects the time scale of the large-scale, non-radial 
pulsations exhibited by the post-shock region (~ 10 
to ~ 20 msec; FigJ^J, while the short-time variations 
(~ msec) of the wave amplitude of the mass flow are 
caused by the convective activity in the hot bubble re- 
gion. The overall rise of the wave amplitude of the neu- 
trinos is a consequence of the steadily increasing rota- 
tion rate (accretion of matter with higher angular mo- 
mentum) and thus growing rotational flattening of the 
nascent neutron star. In the non-rotating model (Fig.[3J) 
and the proto-neutron star model (Fig.0J) the temporal 
variations of the wave amplitudes of the neutrinos are 
caused by large-scale, localized neutrino-emitting down- 
flows through the neutrino-heated bubble, and by global 
changes in the location and shape of the convective emis- 
sion regions, respectively. 

4.2. Wave spectra and signal-to-noise ratios 
In Figs.EltoElwe show the (logarithm) of the quan- 
tity ^(z^)!;/ 1 / 2 appearing in the numerator of Eq. 
Each figure further gives /i rm s(^), the rms noise level 
(see denominator of Eq.EJl of a LIGO I and LIGO II sin- 
gle interferometer (cf. Sect. 3) as a function of the fre- 
quency v, the LIGO II curve being the solid curve of 
Fig-H Both quantities determine the signal-to-noise ra- 
tio (S/N), which is given for all models and both an- 
tennae in Tabled using again the solid sensitivity curve 
of Fig.^ We remind the reader that the S/N values are 
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Table 1. Signal-to-noise ratios' 1 for LIGOI and II. 



Fig. 7. — Spectral energy distribution of the quadrupole radiation 
due to convective mass flow and anisotropic neutrino emission for 
the rotating delayed explosion model sl5r of Buras et.al 12003) as 
a function of frequency of the emitted gravitational radiation for a 
supernova at a distance of 10 kpc. 
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Fig. 8. — Same as Fig.|2]but showing only the contribution due 
to anisotropic neutrino emission. 
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Fig. 9. — Same as Fig.[7|but for the bounce signal only, and a 
supernova at a distance of 3 kpc. 



computed for the complete LIGO I and LIGO II detectors 
as explained in Sect. 3. 

The spectrum of the complete wave train of model sl5r 
caused by both anisotropic mass flow and neutrino emis- 
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sion is displayed in Fig. for a supernova at a distance 
of 10 kpc. It rises to a broad peak between ~ 600 Hz and 
~ 700 Hz and then declines rapidly at higher frequen- 
cies. The figure indicates that while the signal-to-noise 
ratio is probably too small for this event to be detectable 
by LIGO I, it could be very well detected by LIGO II. 
This is confirmed by the signal-to-noise ratios computed 
from Eq. (J2J (using the solid sensitivity curve of Fig.0 , 
which are 3.7 and 67 for the LIGOI and LIGO II de- 
tectors, respectively (Table©. Even at a distance of 
100 kpc the signal-to-noise ratio is still sa 7 for LIGO II, 
it., it is probably sufficiently large for a detection of 
the event (Flanagan & Hughes 1998). Comparing Fig.0 
with Fig. [S] which shows the spectrum resulting from 
anisotropic neutrino emission alone, one recognizes that 
the low-frequency part of the spectrum (below ~ 100 Hz) 
is dominated by the contribution of the neutrinos, be- 
cause their quadrupole wave amplitude varies on much 
longer time scales than that of the mass flow (Fig.[2j) . 

When we Fourier transform just the bounce signal (in- 
sert in Fig.|21), ie., the wave signal up to 15msec past 
core bounce, we find that its spectrum peaks at a some- 
what higher frequency of ~ 900 Hz (Fig.^J) than that of 
the complete signal. The signal-to-noise ratios drop to 
0.7 (LIGOI) and 14 (LIGO II) even when the source is 
located at a smaller distance of 3 kpc (Table© • 

The spectrum of the non-rotating model sllnrl80 
fFig. H0|) is qualitatively quite similar to that of the ro- 
tating model sl5r (Fig.©, but the corresponding signal- 
to-noise ratios are much smaller for a source located at 
the same distance (10 kpc): S/N = 0.6 (LIGOI) and 
S/N = 13 (LIGO II), respectively (Table©. Thus, while 
a non-rotating supernova (model sllnrl80) is detectable 
with LIGO II up to a distance of ~ 15 kpc, a massive 
star rotating at a rate in the ballpark of predictions by 
state-of-the-art stellar evolution models will produce a 
gravitational wave supernova signal which can clearly be 
seen by LIGO II throughout the Galaxy and in the Mag- 
ellanic Clouds, and by LIGOI if the supernova occurs 
less than ~ 5 kpc away from Earth, a distance limit that 
encompasses all historical Galactic supernovae. As in 
case of the rotating model, the low frequency part of the 
spectrum of model sllnrl80 is dominated by the contri- 
bution from anisotropic neutrino emission (Fig.llljl be- 
cause of the slower temporal variation of the correspond- 
ing quadrupole wave amplitude (Fig.©. 

According to our results the "long time" (> 1 sec) grav- 
itational wave signal produced by all core collapse super- 
nova events due to convective activity inside the proto- 
neutron star is detectable with LIGO II (S/N = 8) up to 
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Fig. 11. — Same as Fi gjure|8 lbut for the delayed explosion model 
sllnrl80 of lBuras etoll 12003V 
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Fig. 12 . — Same a s Fig.Qbut for the proto- neutron star model 
pnsl80 of lKeTil <1997h . 



a distance of lOkpc fFig.lT2l and Tabled . The spec- 
trum displays a narrow dominant peak at ~ 1.3 kHz, 
and two additional, weaker broad maxima centered at 
1 kHz and ~ 350 Hz. Anisotropic neutrino emission again 
dominates the spectrum at frequencies below ~ 100 Hz 
(Fig. ll3|> . We point out here that the convective activ- 



ity in the proto-neutron star may continue beyond the 
time period covered by our simulation, perhaps lasting 
more than 10 sec until the deleptonization is complete 
and the convective motion gets damped by viscous ef- 
fects. Assuming that the amplitude and the spectrum of 
the gravitational wave emission does not change consid- 
erably during this time, the signal-to-noise ratios quoted 
for this model in Tabled may increase by a factor of 
> 10, as they depend linearly on the size of the Fourier 
amplitude which grows proportional to the length of the 
signal. 

The signal-to- noise ratios given in Tabled f° r the com- 
plete tuned Advanced LIGO depend somewhat, but not 
extremely, on the choice of the sensitivity curve. In 
case of model sl5r the signal-to-noise ratio increases only 
slightly from 67 to 71 and 70 when using the dotted and 
dashed curves of Fig.^ respectively. The corresponding 
signal-to-noise ratios for model sllnrl80 are 13, 17 and 
10, and those for model pnsl80 are 8, 5 and 8. Hence, 
for the latter two models tuning of the senstivity curve 
for source optimization can increase the S/N values by 
up to 70%. 

4.3. Amount of Radiated Energy 

The energy radiated in the form of gravitational 
waves is given by (see eg., Eqs. [22], [26] and [28] in 
IZwerger fc Mulled 1)1997(1 ) 

c 3 1 f°° fdA^(t)\ 2 

4ikf^r-^, (5) 

where A§f is the Fourier transform of Af 2 , and c and G 
are the speed of light in vacuum and Newton's gravita- 
tional constant, respectively. We used both expressions 
above to compute the radiated energy, and found for all 
three models very good agreement with numerical differ- 
ences only in the few percent range. 

Our analysis implies that non-spherical core collapse 
supernovae rotating with rates in the ballpark of pre- 
dictions by state-of-the-art evolutionary calculations of 
massive stars emit an energy of a few lO _9 M c 2 (Mq 



8 



Miiller et.al 



is the solar mass) in the form of gravitational radia- 
tion (Tabled . This number might increase for stars 
more massive than the considered 15 Mq progenitor, 
and is reduced by roughly one order of magnitude when 
non-rotating cores or the emission of convective proto- 
neutron stars are considered. Anisotropic neutrino emis- 
sion does not contribute much to the amount of radi- 
ated energy, as the corresponding quadrupole wave am- 
plitudes show a slower and smaller temporal variation 
than those due to the non-spherical mass flow. 

5. DISCUSSION AND CONCLUSIONS 

We have presented the gravitational wave signature of 
core collapse supernovae using the presently most real- 
istic input physics available. Our analysis is based on 
state-of-the-art progenitor models of rotating and non- 
rotating massive stars, whose axisymmctric core collapse 
is simulated by integrating the hydrodynamic equations 
together with the Boltzmann equation for the neutrino 
transport including an elaborate description of neutrino 
interactions, and a realistic equation of state. We have 
computed the quadrupole wave amplitudes, the Fourier 
wave spectra, the amount of energy radiated in form of 
gravitational waves, and the signal-to-noise ratios for the 
LIGOI and LIGOII detectors resulting both from non- 
radial mass motion and anisotropic neutrino emission. 

The results show that the dominant contribution to 
the gravitational wave signal is not produced by the core 
bounce itself, but instead by the neutrino-driven convec- 
tion in the post-shock region lasting for up to several hun- 
dred milliseconds post bounce. This finding is qualita- 
tively different from that of previous core collapse simula- 
tions which have been performed with greatly simplified 
parameterized models involving (i) a polytropic equation 
of state, and/or (ii) a simplified description of weak inter- 
actions and neutrino transport or none at all, and/or (iii) 
parameterized pre-collapse stellar models, and/or (iv) 
Newtonian gravity. We further find that the contribution 
to the signal from anisotropic neutrino emission is dom- 
inant at frequencies below ~ 100 Hz, while non-spherical 
mass flow dominantly contributes in the range between 
~ 300 Hz and ~ 1200 Hz. Besides the neutrino-driven 
convection in the post-shock region, the Ledoux convec- 
tion inside the deleptonizing proto-neutron star also gives 
rise to an important gravitational wave signal from core 
collapse supernovae. 

For a stellar iron core rotating with a central angu- 
lar frequency of 0.5 s" 1 , which is relatively small com- 
pared to the rotation rates typically assumed in pa- 
rameter studies b a sed on rotating polytropes (see, eg., 
IZwerger fc Mullen 119971) : IDimmelmeier. Font fc Mullerl 
l|2002fl : IOtt etalll|2003|) T. but which is of the order of the 
largest rotation rates predicted by sta te-of-the-art evolu - 
tionary calculations of massive stars ijHeger e£.all l2003) , 
the core bounce signal would be detectable (S/N > 7) 
with LIGO II for a supernova up to a distance of ~ 5 kpc. 
However, the signal caused by neutrino-driven convection 
is observable with LIGO II up to a distance of ~ 100 kpc, 
and even with LIGO I to a distance of ~ 5 kpc. 

Due to the smallness of the assumed initial rota- 
tion rate we do not expect any significant growth of 
non- axisymmctric instabilities (eg., bar modes) even af- 
ter core collapse when angular momentum conservation 
leads to a spin-up of the core. For the same reason 



there also occurs no centrifugal hang-up, but a regular 
bounce at about nuclear matter density. Consequently, 
the gravitational wave (bounce) signal is not of type II 
consisting of several distinct spikes, but of type I with 
one prominent spike at b ounce followed by a ring-down 
HMorichmever f±r;ill1 991ft. 

As far as the initial state of rotation is concerned, 
our model represents not an extreme, rare case, but 
is in accord with more typical iron core conditions as 
expected from the evolution of rotating massive stars. 
We point out here that most models used in previ- 
ous investigations have too large rotational rates in 
co mparison to those observed for n ewborn pulsars (see 
eg., iKotake. Yamada fe Satol p00 3U. Even our initially 
quite slowly rotating model will produce a pulsar spin- 
ning with a period, P ro t, between 1 and 2 msec (depend- 
ing on its final radius of 10-15 km and assuming angular 
momentum conservation in the neutron star after the end 
of our simulations), which is not far from the breakup 
limit. Millisecond pulsars have rotational energies 
- 51O 52 erg(P rot /lmsec)- 2 (i?/lOkm) 2 M/(1.5M ). If 
such rapid rotators have to be slowed down to match the 
rotation rates of observed young puls ars, which is diffi- 
cult to achieve on short time scales l|Wooslev fc Hegerl 
2003), the question arises where all this rotational en- 
ergy has gone, if it is not observable in kinetic energy or 
electromagnetic radiation from the supernova or pulsar? 
R-modes do also n o longer appear as a likely energy drain 
■l\rms rl. gin-Mil). 

If the core is non-rotating (model sllmT80) its grav- 
itational wave signal due to aspherical motion in the 
neutrino-heated bubble and due to anisotropic neutrino 
emission is detectable (S/N > 8) with LIGOII up to a 
distance of ~ 15 kpc, while the signal from the Ledoux 
convection (model pnsl80) in the deleptonizing proto- 
neutron star can be measured at least to a distance of 
~ 10 kpc. Both kinds of signals are generically produced 
in any core collapse supernova. 

As our models are axisymmetric, the results may 
change when 3D effects associated with the convection 
inside the proto-neutron star and in the hot bubble 
region arc included. Both phenomena are genuinely 
three dimensional. In this respect the simulations of 
iFrver fc WarrerJ l|2002jh and their analysi s of the grav- 
itational wave emission ijFrver et.a I 2003) may provide 
an answer to this uncertainty in our results. These au- 
thors find that the size and dynamics of the convective 
structures are similar to those in axisymmetric simula- 
tions. 

Finally, we want to mention another mechanism which 
causes a gravitational wave sig nal in core collapse ex- 
plosions. Recent simulations bv lScheck et.al\ 1)2003(1 in- 
dicate that the proto-neutron star is accelerated during 
the explosion at a rate of up to ~ 3000 km/sec 2 , result- 
ing in large kick velocities in some models. The acceler- 
ation can be caused either by asymmetric accretion onto 
the proto-neutron star, or by the pull of the gravita- 
tional potential of matter distributed asymmetrically in 
the surrounding hot bubble region. Due to the erratic 
behavior of the accretion downflows the acceleration can 
temporarily change sign leading to an additional oscilla- 
tory component superimposed on or preceding the steady 
acceleration component. 

An order of magnitude estimate of the resulting scale 
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of the gravitational wave signal can be obtained by ap- 
proximating the proto-neutron star by a point mass M 
(initially at rest), which eg., moves in positive z-direction 
starting at z = 0. If the point mass experiences a 
constant acceleration b for a time t, its displacement is 
given by z(t) — 6/2 t 2 . The only non-vanishing (Carte- 
sian) component of the time-dependent quadrupole mo- 
ment (measured with respect to the origin of the coor- 
dinate system) of such a moving point mass is Q zz = 
2/3 Mz 2 (t) = l/6Mi)¥. The lowest (quadrupole) or- 
der contribution to the gravitational wave luminosity 
L = d£cw /dt of such a moving point mass is du e to the 
temp oral variation of its quadrupole moment ijThornel 
1980), and is given by L kick = G/(5c 5 ) ((d 3 Q zz /dt 3 ) 2 ) = 
16G/(5c 5 ) M 2 bH 2 , which can be written as 



L kick = 2.01(T 19 ( -— 



M 

Me 



b\t 2 



(6) 



where 62 is the acceleration in units of 10 2 km/sec 2 , 
and t the time of acceleration, measured in seconds. 
If the motion of the point mass (initially at rest at 
z = 0) is caused by a periodic acceleration of ampli- 
tude A and time scale t, ie., d 2 z(t)/dt 2 = A cos(27rf/r), 
the resulting gravitational wave luminosity is given by 



L osc = 2G/(9tt 2 c 5 ) M 2 A 4 t 2 , or by 





[M c 2 l 




sec 



(7) 



where A3 is the acceleration in units of 10 3 km/sec 2 , 
and T—i the oscillation time scale, measured in units 
of 10 _1 sec. From their simulations IScheck et.al\ (2003) 
infer b% < 5, t & 1, A3 < 3, and r_i ~ 1 for the 
proto-neutron star of mass M ~ Mq. Hence, the 
gravitational wave signal from the acceleration giving 



rise to the kick of the proto-neutron star is quite small 
(^kick < 10~ 16 [MqC 2 /sec]), and the signal from the os- 
cillatory motion of the proto-neutron star is even smaller 
(^osc ;$ 10~ 17 [M Q c 2 /sec]). The former signal is ex- 
pected to have a frequency of typically ~ Hz, and the 
latter one should be dominant at about 10 to 20 Hz. 

Both signals are very small compared to the sig- 
nal caused by the bounce of a rotating core, because 
the accelerations encountered during core bounce are 
much larger. The accelerations, which cause the varying 
quadrupolar deformation of the rotating core, are of the 
order 10 6 km/sec 2 , ie., they are roughly a factor of 10 3 
larger than those producing the kick or the oscillatory 
motion of the proto-neutron star. These accelerations 
act for only ~ 10 msec (the damping time scale of the 
post-bounce ringing; see eg., the insert in Fig.[2Jl, which 
is a factor ~ 100 (~ 10) shorter than t (t_i), and involve 
a mass of only ~ 0.1 M Q , which is a factor ~ 10 smaller 
than M. Nevertheless, as L scales bi-quadratically with 
the acceleration and only quadratically with the mass 
and time scale (see Eqs.[£j]and[3), the gravitational wave 
luminosity of rotational core bounce exceeds that of the 
proto-neutron star (kick and oscillatory) motion by a fac- 
tor ~ 10 6 to - 10 8 . 
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